An alternative non-Markovianity measure by divisibility of dynamical map 
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Identifying non-Markovianity with non-divisibility, we propose a measure for non-Markovinity of 
quantum process. Three examples are presented to illustrate the non-Markovianity, measure for non- 
Markovianity is calculated and discussed. Comparison with other measures of non-Markovianity is 
made. Our non-Markovianity measure has the merit that no optimization procedure is required and 
it is finite for any quantum process, which greatly enhances the practical relevance of the proposed 
measure. 
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I. INTRODUCTION 

A quantum process is said to be Markovian if the fu- 
ture states of the process depends only on the state of 
present time. In contrast, dependence on past states 
should then be a characteristic feature of non-Markovian 
processes. With the development of technology to ma- 
nipulate quantum system, the quantum non-Markovian 
process has attracted increasing attention in recent years 
[l|-l6j . On one hand the inevitable interaction of a quan- 
tum system with its environment leads to dissipation of 
energy and loss of quantum coherence, on the other hand 
the quantum system may temporarily regain some of the 
previously lost energy and/or information due to non- 
Markovian effects in the dynamics. This motivates the 
study on the non-Markovianity and a measure for the 
degree of non-Markovianity is indeed needed. 

Several approaches are proposed to quantify non- 
Markovianity, including the measure based on the in- 
crease of trace distance [Zi], the measure by quantifying 
the increase of entanglement shared between the system 
and an isolated ancilla and the divisibility of the dynami- 
cal map 8] , the measure based on the decay rate of master 
equation itself 9] , and the measure through the Fisher in- 
formation flowjlflj. Although several approaches to quan- 
tifying the non-Markovianity are proposed, the definition 
of non-Markovianity still remains elusive and, in some 
sense, controversial |llj. 

It has been proven that all divisible dynamical maps 
are Markovian, this divisibility property holds for a larger 
class of quantum processes than those described by the 
Lindblad master equation, for example, the time-local 
master equation with positive decay rates. This indi- 
cates that the divisibility may be a good starting point 
to quantify non-Markovianity. In this paper, we will pro- 
pose a measure for non-Markovianity based on the divisi- 
bility of dynamical maps, three dynamical maps are pre- 
sented and the corresponding non-Markovian measures 
are calculated. These results suggest that the measure 



can capture the feature of non-Markovian dynamics, and 
provide an easy way to calculate the non-Markovianity. 

The paper is organized as follows: In Sec. II, we dis- 
cuss the non-divisibility of dynamical map and the non- 
Markovianity of this map. In Sec. Ill, we introduce a 
measure for non-Markovianity. Three examples are pre- 
sented and discussed in Sec. IV. Sec.V summarizes our 
results. 



II. NON-DIVISIBILITY AND 
NON-MARKOVIANITY 

In quantum Markovian process, the future state of the 
quantum system depends only on the state of present 
time. However, writing this statement in a precise math- 
ematical representation is not an easy task. Instead, we 
use the following description for quantum Markovian pro- 
cess. A quantum evolution is Markovian if it is an ele- 
ment of any one-parameter continuous completely posi- 
tive semigroup [12j . The quantum evolution governed by 
the master equation 



dp ^ 



(1) 



is an example, where £ is a time-independent generator 
of the well-known Lindblad form, 

Cp - -t[H,p] 

+ ^7a(^aPl4^-^K.V„p-ipF„V„) (2) 

a 

with 7c( > 0. This generator leads to completely posi- 
tive trace-preserving maps A(t) — e'^* and it satisfies the 
composition law. 



A(ti-f ta) = A(t2)A(ti). 



(3) 



If a dynamical map can be written in this decomposition 
with both A(t2) and A(ii) being completely positive, the 



2 



dynamical map is called divisible. This composition law 
can be extended to a general case, where the generator 
in Eq. ([2]) is time-dependent, namely. 



dp 
IE 



C[t)p 



with 



C{t)p = -i[H{t),p]+Y.^^{t){Vo,{t)pVj{t)- 
\vJ{t)V^{t)p-]^pVj{t)V^{t)), 



(4) 



(5) 



where ^a{t) > 0. This is known as time-dependent 
markovian [l^. The solution to Eq. (|4]) can be writ- 
ten in terms of the two-parameter family of dynamical 
maps A{t2,ti) {t2 > ti > 0). The composition law corre- 
sponding to Eq. ^ becomes 



A(t2,0) = A(t2,ii)A(ti,0), 



(6) 



and the map A(i2,ii) can be written as A(t2,ii) — 

Te-'*i ^ ' , r is the chronological operator. The com- 
position law (divisibility of the map) implies that the 
dynamical map A(i2,^i) 



A(t2,ii) = re^'''^(*')''*' {t2>h>0), 



(7) 



transforming a state at ti into a state at t2 (for sys- 
tems governed by time independent master equation ([T]), 
A{t2^tl) — K{t2 — ti,Q) — A(<) ) must be trace-preserving 
and completely positive, regardless of which dynamics 
it describes, e.g., it is from the time-dependent master 
equation or the time-independent master equation. Note 
that the starting time ti is not zero. 

A measure for non-Markovianity should quantify the 
deviation of a dynamical map from Markovian evolution. 
Noticing that when a dynamics is non-Markovian, the 
dynamical map A(t2,ii) may not be completely posi- 
tive, we may use the non-divisibility to quantify the non- 
Markovianity. In fact, this is the underlying reason that 
the trace distance can increase 0], and the system gains 
entanglement with an isolated ancilla f§|. 

It is worth stressing that there is no contradiction be- 
tween the requirement on non-completely positivity and 
that on physics. Consider a quantum evolution in a time 
interval (0,^2), we always have A (t2, 0) — A(i2, ti)A(ti, 0) 
due to the time continuity. For A(i2, 0) to be a dynamical 
map, it is required that A(t2, 0) must be completely pos- 
itive, however, A(t2,ii) may not be completely positive. 
Therefore these two-parameter maps in non-Markovian 
dynamics do not generate a quantum dynamical semi- 
group. Then one may wonder: does there exist a A (^2, ii) 
that it is not completely positive but A(t2,0) does? The 
answer is yes. First, a wide range of non-Markovian pro- 
cess can be described by time-local master equations via 
time-convolutionless projection operator |13l - [l8j . Sec- 
ond, it has been shown that any quantum dynamics de- 
scribed by memory kernel master equation may be writ- 
ten in a time-local form [19]. Note that the decay rates in 



these time-local master equation are different from that 
in Eq. ([S]), they can be negative. With this time- local 
master equation, the dynamical map with non-zero start- 
ing time in Eq.([7]) may violate the complete positivity due 
to the negative decay rates. This implies that the non- 
complete positivity of the map A(i2,^i) is an essential 
feature of non-Markovian process. The time-dependent 
decay rate may be negatively infinite at some points of 
time, where the reviving of population or regaining of 
quantum coherence happen [3, [l3l ■ We call these points 
of time singular points tg. When ti = ts or t2 — tg , 
A{t2,ti) may not exist. However, we can use A{t2,ti) in 
the limit that ti — >■ tg instead of A{t2,ts). It is conve- 
nient to discuss A{t2,ti) with a specific time-local master 
equation, but this is not necessary. 



III. MEASURE FOR NON-MARKOVIANITY 

To construct a measure for non-Markovianity, we re- 
sort to the Choi-Jamiolkowski isomorphism [20, HH, it 
asserts that a linear map A -.Md — > Md is isomorphic to 
the Choi matrix. 



(8) 



where \i) are orthogonal bases. An familiar form of the 
Choi matrix is 



PA = {A(E,lM{cf>\, 



(9) 



where \(j>) is the maximally entangled state |0) = 
^ X^iLi K) ® K)i I is an identity map acting on the an- 
cilla, and Pa is the normalized Ca- It turns out that A is 
completely positive if and only if pA is positive semidef- 
inite. In other words, the sufficient and necessary con- 
dition of non-complete positivity is the negativity of pA- 
Hence, the sum of negative eigenvalues of pA can be taken 
as a measure for the non-complete positivity of the dy- 
namical map. However, the summation may sometimes 
be an infinite value in some models due to singular decay 
rates, this suggests to use a normalized quantity 



Ncp = arctan(— Xk 



(10) 



as a measure for non-complete positivity of the map A, 
where is the k ^th negative eigenvalue of pA- Clearly 
if PA > 0, Ncp = 0. Then we have < Ncp < f . Based 
on the complete positivity property of the map A(ti, ^2), 
a measure of non-Markovianity has been proposed [8|. 
This measure is different from ours in that we use the 
an averaged negativity of the map A{ti,t2) to quantify 
the non-Markovianity. Moreover, insightful examples are 
presented to shed light on the non-Markovianity measure. 

To calculate Eq. (|10p with a given time-local mas- 
ter equation, the exact form of A{ti,t2) is not necessary. 
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What we need is to extend the time-local master equa- 
tion from one system to two systems, taking an isolated 
ancilla attached to the system. The Hilbert space is ex- 
tended from Hd to Hd ® Hd accordingly. All operators, 
say O, are replaced by 0®1. By this extension, we get a 
new master equation, which describes the system and the 
ancilla. The system evolves in the same manner as before, 
while the anciila is isolated from both the system and en- 
vironment. The master equation can be solved starting 
from |(/)) at time ti and the state at time t2 {t2 > ti) can 
be obtained. 

We aim at finding a measure NM for non- 
Markovianity which captures the feature of non-complete 
positivity of all possible A{t2,ti). Note that Ncp is a 
function of ii and t2. Let S count the number of Ncp in 
all time interval with Ncp > 0, i.e. 



S 



dti I dt2 c{t2,ti), c - 



1, if Ncp > 0, . 
0, if Ncp = 0} > 



If 5 = 0, i.e., all A{t2,ti) (for any ti and t2, as long 
as t2 > ti) are completely positive, the non-Markovinity 
NM should be defined to be zero. If 5* > 0, we define 



NM 



lim 



J^dhJ^^dt2 Ncp {t2,ti) 
/o^'^^l C{t2,ti) 



(12) 



as a measure of non-markovianity. This can be under- 
stood as an averaged non-complete positivity of all the 
non-completely positive maps in interval t = (0,+oo). 
Therefore, from < Ncp < f , we have < NM < f . 

The upper limit of the integral T in the definition Eq. 
(1121) is taken to be infinite. However, the distribution 
of Ncp on (^1,^2) plan is often periodic or is limited in 
a small area, this suggests that the integration can be 
taken merely in one period of time or taken in a finite 
area. When Ncp is neither periodic nor limited in a small 
region, the upper limit T in Eq. ([T^ should be large 
enough to get a convergent NM. The definition can be 
written into a simple form. 



NM = EiNcpiAN)), 



(13) 



where A^r represents all the non-completely positive 
maps and E means expectation value. Therefore, NM 
can be numerically calculated by averaging a large num- 
ber of non-completely positive maps with equal weight 
(or randomly) in a representative area as discussed above. 
The representative area means one period in time or a 
limited region where Ncp > 0. 

To illustrate the measure of non-Markovianity, we 
present three examples in the next section. We work 
in the interaction picture for simplicity to calculate the 
measure, since unitary transformation does not change 
the eigenvalues of pA as well as Ncp of A{t2,ti), hence the 
non-Markovianity measure under unitary transformation 
remains unchanged. 



IV. EXAMPLES 

A. damping J-C model 

The first example is a two-level system coupling to a 
reservoir at zero temperature. The reservoir consists of 
infinite number of harmonic oscillators that is also re- 
ferred in the literature as the spin-boson model. This 
model is exactly solvable [l3| • The Hamiltonian for such 
a system reads. 



H = Hn 



H, 



with Hq = hjjJt:,aj^a^ -|- hujkb^bk 



Hi = a+B + cr^B^ 



(14) 



, (15) 
, (16) 



where B ~ gkbk- The Rabi frequency of the two-level 
system and the frequency for the k — th harmonic oscilla- 
tor are denoted by ujq and Wfe, respectively. and are 
the creation and annihilation operators oi k — th oscilla- 
tor, which couples to the system with coupling constant 
9k- 

Assuming the system and the reservoir are initially un- 
correlated, we can obtain a time-dependent master equa- 
tion in the interaction picture. 



1 



7(i)(cr pcT+ 
= -2Im[^ 



-a - a p- -pa a 



(17) 



where s(t) = -2Im[|^] and jit) = -2Re[^]. m 

plays the role of Lamb shift and 7(t) is the decay rate. 
Both n{t) and j{t) are time-dependent. c{t) is deter- 
mined by c{t) — ~ Jq f{t — T)c{T)d{T), where f{t — t) — 
J dLuJ{uj)exp{i{ujQ—uj){t — T)) is the environmental corre- 
lation function. In the derivation of the master equation, 
the reservoir is assumed in its vacuum at i = 0. 
Consider the following spectral density. 



J{uj) = -- 



(18) 



where 70 represents the coupling constant between the 
system and reservoir, A defines the spectral width of the 
coupling at the resonance point ujq. For the spectral den- 
sity (UHl), we have s{t) = 0, c{t) = coe-^*/2[cosh(f ) -|- 

with d = 



sinh(f )], and j{t) = j 



27oAsinh(dt/2) 



cosh((it/2)-hAsinh(dt/2) 



\/ — 270A. Note that Eq. ((T7| is derived without any 
approximations, hence it is non-Markovian and it exactly 
describes the dynamics of the open system. 

It is well-known that A characterizes the correlation 
time Tfl of the reservoir through r^j — . The time 
scale Tg on which the state of the system changes is given 
by Ts = 7(7^- So the degree of non-Markovianity should 
be relevant to the rate i? = th/ts- Namely, when R is 
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FIG. 1: Ncp and non-Markovianity in the damping J-C 
model, (a) Ncp versus ti and At (in units of 1/A) at i? = 5. 
Three oscillations are shown, (b) NM as a function of R. 
Clearly, NM monotonically increases with R. 



very small, the evolution is Markovian, when tr is com- 
parable with Ts , the memory effect of reservoir should be 
taken into account, the dynamics of the open system is 
then non-Markovian. 

Let us first analyze the non-Markovianity of the dy- 
namics by examining 7(t). For R < i, 7(t) is always 
positive and it is a monotonically increasing function of 
time, all A{t2,ti) are completely positive, hence the dy- 
namics is Markovian. When i? > i, 7(t) is a periodic 
function of time, it takes negative values sometimes. In 
particular, 7(t) has discrete singular points where the 
upper level gains population, this is a typical feature of 
non-Markovianity. 

Now we see if our measure can capture all these fea- 
tures of non-Markonianity. In order to apply our mea- 
sure, we have to extend the time-local master equation 
to the compound system, i.e., the operators a± in Eq. 
p?| is replaced hy a±(S)l with I being the ancilla's 2x2 
identity operator. For a given time interval (ii,t2), a 
straightforward calculation yields 



Ncp{ti,t2) 



j arctan(i(^-l)), 



0, 



> 1 

< 1, 

(19) 

where c{t) was defined below Eq. (|17p and t2 > ti > 0. 

For a typical non- markovian case (R — 5), we plot Ncp 
as a function of ti and At in FiglIJa) (At = t2—ti, since 
t2 > ti, we use At instead of t2 for convenience). This 



FIG. 2: Ncp and non-Markovianity in the detuning J-C 
model, (a) Counter plot of Ncp as a function of ti and At 
(in units of 1/A), with A = 10 and 70 — 0.3 (in units of A), 
(b) NAI as a function of A. Other parameters chosen are 
the same as in (a) . The maximum non-Markianity arrives at 
about A = 6. 



plot shows the non-zero area and its value of Ncp versus 
ti and ti — t2. As Ncp is a periodic function of ti, and 
the area where Ncp > decays very fast with At. NAI 
can be given by averaging all Ncp in one period, it yields 
NM = 0.835. 

As expected, when R < 0.5, NM = 0, and when R > 
0.5 the non-markovianity is finite. We plot the measure 
of non-Markovianity with different R in FigHJb). Here 
R is chosen from 0.5 to 10. Note that when R = 0.5, 
7(t) does not exist due to the zero deirominator. Our 
result for R = 0.5 is obtained at R — 0.5 + e with e 
an infinitesimal positive number. Intuitively, the larger 
R is, the stronger the non-Markovianity. The results in 
FigHJb) demonstrate that this is indeed the case. 



B. J-C model with detuning 

The second example is similar to the first one, but 
here we consider the system in a cavity whose center 
frequency is detuned from the system Rabi frequency wq. 
The dynamics in the interaction picture is governed by 
Eq. PT)) . but s{t) and j{t) are determined by the Lorenz 
spectral density 



J(w) 



1 



70AV2 



27r {ujQ 



;)2 + A2' 



(20) 
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where A denotes the detuning. 

We extend Eq. PT|) to a compound system by intro- 
ducing an ancilla as we did in the first example, then 
we calculate Ncp{ti,t2) numerically. We plot Ncp as 
a function of ti and At in Fig[2][a) for a typical non- 
Markovian case. We see that Ncp{ti,t2) is finite in 
contrast to infinite values in the same region for the 
first example. On the other hand, Ncp of all A(i2,ti) 
(with different ti and are far less than ^, indicating 
weaker non-Markovianity in comparison with the first ex- 
ample. Finally, from Eq. (HH), or Eq. we have 
NM = 3.66 X 10"'' in this case. 

Now we discuss the dependence of non-Markovianity 
on the detuning. We plot NM in FiglJJb) with different 
A. We find that NM appears non-zero at about A = 
4, it first increases then decreases with A. This result 
is similar to that in Q, where the non-Markovianity is 
measured by the decreases of trace distance. 



C. a two-level system coupling to a finite spin bath 

In the third example, we consider a central spin-^ cou- 
pling to a bath of N spin-|s. The interaction Hamilto- 
nian is, 



N 



(21) 



fe=i 



Pit) 



Pii 



P21C0& (^) 



Pl2 COS 

P22 



N(2At\ 



(22) 



In terms of dynamical map, the dynamics can be rep- 
resented as, A{t,0)p = i(l - cos^{^j)a^p(T^ + i(l + 

cos''^(-^^))p. This is equivalent to the following master 
equation. 



(23) 



where C{p) — <JzP<Jz — p, and the time-dependent ^{t) ~ 
^ tan(^^). This example is discussed in several papers 
as a classical example to quantify non-Markovianity, and 
the non-Markovianity is infinite^, i8|]. By our definition 
of non-Markovianity, it is finite. This allows us to estab- 
lish a relation between non-Markovianity NM and the 
number of spin N in the reservoir. 

It is easy to calculate pA (defined in Eq.®) by Eq. 



(A(t2,ti)®l)|V^)(V'l - 



0.5 0.5fc' 





0.5fc 0.5 



(24) 
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where Ak — A/yN represents the coupling constants. 
The non-Markovianity of the central spin in this model 
is discussed in [ij] . Assume the initial state of the whole 
system is Ps(0) ® (jnI): i-e., all spins in the reservoir with k 
are in a maximal mixed state. The density matrix of the 
central spin at time t takes. 



FIG. 3: Ncp and non-Markovianity in the N-spin bath model, 
(a) Ncp as a function of ti and At for A'^ = 1. t is in units 
of It shows that Ncp is a periodic function of ii and t2, 
two periods are shown here, (b) NM with different N. Non- 
Markovianity increases with N and tends to j when the limit 
N — !> +00 is taken after T — !> oo. 



JV / 2At2 ^ 



(^2 > ti > 0), then we have 



Ncpituts) 



arctan(i( 
0, 



1)), 



(|fc| > 1) 
-1 < fc < 1). 



(25) 

The non-completely positivity of the map A(t2, ti) can 
be examined by plotting Ncp in (ti,t2) plan. FiglSja) 
shows this result taking only one spin as the environment 
(for > 1, the results is similar). We find that Ncp is a 
periodic function of both ii and Ai. Whenti = (2n-|-l)^ 
and At {n — 0,1,2---), Ncp — ^ reaches its 

maximum. The non-Markovianity measure in this case 
is NM = 0.505. 

Now we try to find the relation between non- 
Markivianity and the number of spins A^ in the envi- 
ronment, the result is plotted in Fig|31^b). We find that 
NM increases with the spin number N . For a very 
large A'^, NM arrives at its maximum. However, the 
coupling strength Ak ~ A/^/N is very small in large 
A'^ limit, this implies that the characteristic time of the 
central spin tends to be infinitely long in this limit, i.e., 
non-Markovian dynamics can be observed only on a long 
time scale. If we are interested in a limited time interval 
(0,t), all dynamical maps A(ti,t2) are completely posi- 
tive, there is no non-Markovian effect. This can be under- 
stood as follows. With a fixed T, one always can choose 
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N that is close to zero. The off-diagonal element 

of the density matrix Eg.lfM]) then takes pi2e~^'* (t2-*i)^ 
indicating that the density matrix in Eq. ([24| describes a 
typical Markovian process. 



optimization is required. Three examples are illustrated, 
which show that the measure can nicely manifest the non- 
Markovianity. We also compare our measure with others 
in the literature and find that it almost agrees with the 
trace-distance based measure. 



V. SUMMARY 

We have presented a measure for non-markovianity 
based on the divisibility of dynamical map. This mea- 
sure has the advantage that it is easy to calculated and no 
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